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We perform a numerical study of the excess of non-abelian gauge invariant gluonic action around 
thermal abelian monopoles which populate the deconfined phase of Yang-Mills theories. Our results 
show that the excess of magnetic action is close to that of the electric one, so that thermal abelian 
monopoles may be associated with physical objects carrying both electric and magnetic charge, 
i.e. dyons. Thus, the quark gluon plasma is likely to be populated by selfdual dyons, which may 
manifest themselves in the heavy-ion collisions via the chiral magnetic effect. Thermodynamically, 
thermal monopoles provide a negative contribution to the pressure of the system. 
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I. INTRODUCTION 

Interpretation of heavy ion experiments at RHIC and 
first-principle lattice simulations suggest that the quark- 
gluon plasma has quite unusual properties Contrary 
to general expectations, at temperatures just above the 
critical temperature, Tc, the plasma looks like an ideal 
fluid rather than a system of weakly interacting quark- 
gluon gas. The reason of this behavior lies in the strongly 
interacting nature of the quark-giuon plasma. 

The magnetic component of the gluon fields may be a 
key player in the unusual dynamics of the plasma 0, H, 
0, H, H, H H, Q . The basic idea underlying the magnetic 
picture can be motivated in the following way As 
it was suggested long time ago in Ref. at low tem- 
peratures the magnetic degrees of freedom exist in the 
form of a (gluo)magnetic condensate. Magnetic charge 
condensation is believed to be responsible for color con- 
finement and this is usually known as the dual super- 
conductor scenario: the existence of such condensation 
and its disappearance at the critical temperature has 
been verified in various lattice QCD studies (see e.g. 
Refs. [n!, [H, [H, [13 or Ref. [H for a review). As the 
temperature increases, the condensate starts to evapo- 
rate and it is destroyed completely at Tc, while the mag- 
netic color degrees are released as thermal states into the 
gluon plasma. The central point is that simultaneously 
the magnetic degrees of freedom are to become physical 
and constitute a real, thermally excited component of the 
gluon plasma 

The magnetic monopoles in the Yang-Mills plasma are 
particlelike objects made of gluons. The monopoles are 
suggested to affect both the transport properties and the 
thermodynamical features of the plasma 0, iS, i9] . This 
suggestion can in principle be checked by numerical sim- 
ulations of Yang-Mills theory on the lattice. The lat- 
tice method provides us with thermalized ensembles of 
monopole trajectories in Euclidean space time. However, 
in order to calculate the contribution of thermal (i.e.. 



real) monopoles to any physical quantity one should be 
able to separate the thermal monopoles from numerous 
virtual monopole loops. Virtual monopoles - which are 
responsible essentially for ultraviolet physics - are known 
to populate the thermal ensembles in numerous quanti- 
ties, so that the separation of monopoles into real and 
virtual parts seems to be a difficult task. 

It was shown in Ref. that thermal monopoles 
in Minkowski space-time are associated with Euclidean 
monopole trajectories that have a nontrivial wrapping 
number on the short (thermal) direction of the Euclidean 
space. Moreover, the density of the monopoles in the 
Minkowski space-time is given by the thermal average 
of the absolute value of the monopole wrapping number. 
The properties of the wrapped trajectories were studied 
numerically using lattice simulations in Ref. I, [3. In 
Ref. [4,] the average wrapping number of the monopole 
trajectories was thoroughly calculated. It was shown that 
this quantity is free from ultraviolet artifacts as it should 
be the case for a purely thermal quantity. The monopole 
density grows with the temperature T. However, the 
parametrical T-dependence of this quantity differs from 
the behavior expected in the case of an ideal gas of the 
monopoles. This is one of the first lattice evidences on 
the nontrivial dynamics of the magnetic monopoles in the 
high temperature gauge theory [J|. 

Monopoles do not appear alone in the plasma. Numer- 
ical simulations have revealed the fact that monopoles are 
a part of chainlike structures [lB| ■ In these structures the 
(anti) monopoles are connected to each other by magnetic 
vortices. The existence of the vortex defects is related to 
the nontriviality of the center of the gauge group. Both 
the monopoles and vortices contribute to the equation of 
state of Yang-Mills theory at finite temperature 0, 0] (a 
short review of the subject can be found in Ref. [8]). 

A calculation of the effect of the vortices on the pres- 
sure of the Yang-Mills plasma can be formulated in 
a rather straightforward and clear way. According to 
Ref. ^ the vortices provide a negative contribution to 
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the pressure in agreement with theoretical expectations 
of Ref. As for the monopoles, the prehminary cal- 
culation of Ref. 0] indicated that they provide a large 
positive contribution. However, the calculation was done 
using a nonlocal method which did not discriminate be- 
tween the real and virtual monopole loops. Below we 
show that the contribution of the thermal monopoles to 
the pressure of gluons is negative as well. 

Thermal monopoles detected by lattice QCD simula- 
tions are of abelian nature. An abelian projection is 
needed to define magnetic monopole currents (see the 
Appendix for details) and the Maximal Abelian Gauge 
(MAG) is usually chosen to do that. The fact that 
both the density and the spatial correlations of thermal 
monopoles measured in the MAG projection are free of 
UV artifacts Q (i.e. they correctly scale to the contin- 
uum limit) points to an underlying physical meaning of 
those objects. Nevertheless, the dependence on the cho- 
sen abelian projection still poses a problem. In order to 
further investigate this issue and better clarify the phys- 
ical properties of thermal monopoles, we have decided 
to measure the correlation of the unprojected gauge in- 
variant nonabelian action density with the positions of 
abelian thermal monopoles. As we shall clarify in the 
following, this is also the essential quantity entering the 
determination of thermodynamic properties of thermal 
monopoles (e.g. their contribution to the pressure). 

We made refined measurements which treat the mag- 
netic and electric parts of the gluon action separately. 
This allowed us to find - as an unexpected byproduct 
- that the excesses of the magnetic and electric action 
densities have almost the same value near the monopoles 
suggesting that the electric and magnetic components of 
these objects have the same significance. We conclude 
that the objects which we identified simply as thermal 
abelian monopoles after abelian projection, may corre- 
spond to more interesting physical objects, actually ther- 
mal dyons, with equal electric and magnetic charges 

The dyons in zero-temperature Euclidean Yang-Mills 
theory were discussed in the literature for a long time 
both in continuum [l^ and lattice [l^, [13] formulations. 
In Ref. [20] smooth dyonic solutions were observed in the 
deconfinement phase of lattice Yang-Mills theory. Here 
we discuss another lattice evidence in favor of existence 
of dyons as real objects in quark gluon plasma. It is 
necessary to mention that these dyons would have a finite 
density which is independent of the ultraviolet cutoff. 

The importance of non-Abelian dyons in the quark 
gluon plasma is difficult to overestimate since these ob- 
jects may play a key role in the so-called chiral magnetic 
effect. This effect - which reveals certain CP-odd cor- 
relations in the quark-gluon plasma ~ was suggested the- 
oretically in Ref. [21]. The chiral magnetic effect leads 
to generation of an electric current along the direction of 
a magnetic field in a nontrivial topological backgrounds 
of gluons. The strong magnetic field is naturally created 
in heavy-ion collisions, and the appearance of the elec- 
tric current is detected in a form of a certain geometric 



asymmetry of electrically charged particles produced in 
the collisions. There are preliminary indications that the 
chiral magnetic effect has been indeed observed by the 
STAR Collaboration in experiments at RHIC [11]. The 
signatures of this phenomenon were also found in recent 
numerical lattice simulations [23j, and discussed in non- 
perturbative holographic approaches of Ref. [23] . In case 
the dyons are really associated to the thermal monopoles, 
as suggested by our results, they may serve as fluctua- 
tions of the topological charge in the quark gluon plasma 
phase. 

The structure of this paper is as follows. In Sec. |TT] 
we review basic thermodynamical relations of Yang-Mills 
theory. In Sec. IIIII we formulate how to estimate the in- 
fluence of the monopoles on thermodynamics and how 
to calculate the electric and magnetic action densities 
around the monopoles. In Sec. IIVI we present our nu- 
merical results. Our conclusion is summarized in the last 
Section, and technical details are given in the Appendix. 



II. THERMODYNAMICS OF GLUONS 



The partition function of the Yang-Mills theory is 



Z= /DAcxp^-^TrG^, 



(1) 



where G^,y — G^^j^f^ is the field strength tensor of the 
non-Abelian field A, and t"" are the SU{N) generators 
normalized in the standard way, Trf^t^ — ^S'^''. 

For a sufficiently large and homogeneous gluon system 
residing in thermodynamical equilibrium the pressure p 
and the energy density e are given by the derivatives 
of the partition function ([1]) with respect to the spatial 
volume, V, and the temperature T, respectively. 



T dlogZ{T,V) _ T 

V a logy ~ V 

T d log Z(T, V) 

V 9iogr ' 



log Z{T,V), 



(2) 
(3) 



Following a standard approach (see, e.g., Ref. [25|), one 
can relate the energy density e and the pressure p can via 
the quantum average the trace of the energy-momentum 
tensor T^i,, 



e{T) = {Tp^e-3p, 
where in Yang-Mills theory 



T., = 2Tr 



G p(jGi/(j fii^G (J pG (J p 



(4) 



(5) 



At the classical level trace of this quantity is zero be- 
cause the bare Yang-Mills theory is a conformal theory. 
However, at the quantum level the conformal invariance 
is broken, and, consequently, the energy-momentum ten- 
sor exhibits a trace anomaly: the quantum average of the 
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trace of the energy momentum tensor ^ is nonzero. The 
trace anomaly is intimately related to the gluon conden- 
sate which breaks the scale invariance of the theory 



(6) 



where /3{g) the Gell-Mann-Low /9-function. We use the 
notation (3 instead of the conventional (3 in order to avoid 
a confusion with the lattice coupling constant 



f3- 



2N 



(7) 



The trace anomaly (jl]) is an important thermodynamic 
quantity because it allows us to reconstruct both the 
pressure and the energy density as follows: 



A f'^ dTi 0(Ti) 
Jo J-i ^1 

Jo ^1 J-i 



(8) 
(9) 



On the lattice the dynamical fields are the gluonic 
SU{N) matrices JXr^ which are identified at the lattice 
links I = {x,^}. The link gauge fields are related to the 
continuum gauge fields A^{a) in the limit of vanishing 
lattice spacing, a 0: 



x-\-afi 



Uxfi = expj ig 



(10) 



l + iagAi_^{x) + O(a^) . 



The lattice spacing a is a function of the lattice cou- 
pling 0. 

The lattice analogue of the partition function ^ is 



Z{T,V)= J i?C/cxp{-/3^5p[[/]} 



(11) 



where the plaquette action density S'p[C/] is usually writ- 
ten in the Wilson form, 

Sp[U] = 1 - ^Re TvUp, Up= UUi, (12) 



ledP 



The spatial volume V — (Lga)^ and the temperature 

1 



T = 



Lta 



(13) 



of the gluonic system are defined by the asymmetric ge- 
ometry of the Euchdean lattice, L^Lt- The shorter direc- 
tion, Lt with Lt ^ Ls, is associated with the imaginary 
time, or, "temperature" direction. The trace anomaly ([6]) 
on the lattice is written as follows 



where {Sp)rp and {Sp)^ are the action densities at T > 
and at T = 0. The subtraction of the T = value 
is needed to remove the effect of quantum fluctuations, 
which lead to ultraviolet (UV) divergency of the quantum 
expectation value. The regularized trace anomaly (|14[) is 
an UV-finite quantity. 

We used the Monte Carlo simulations on the lattice 
to evaluate the influence of the monopoles on the trace 
anomaly, and, consequently, on the equation of state in 
the thermal Yang-Mills theory. We describe setup of our 
numerical simulations in Appendix [XI 



III. LOCAL ACTION DENSITIES 

The (chromo) magnetic monopoles in the Yang-Mills 
plasma are particlelike objects made of gluons which are 
generally believed to be responsible for color confinement 
in QCD. The monopole confinement scenario is reviewed 
in Ref. [l^. On the lattice the monopoles appear in the 
form of the closed trajectories. The monopole loops may 
be either wrapped around the short (imaginary time) di- 
rection of the lattice, or they may form shrinkable loops 
with zero winding number. The monopole configura- 
tions with nontrivial winding are associated with thermal 
monopoles while the other trajectories correspond to 
virtual loops. The properties of the winding trajectories 
were studied long ago in Ref. [Hill]. A recent study Q 
has revealed that the thermal monopole density. 




[mo{x,l 



(15) 



identified as a density of the average wrapping number 
A^wrap['7io(^, t)], shows a very good scaling towards a 
continuum limit. This result confirms the physical nature 
of thermal, or real, magnetic monopoles in the Yang-Mills 
plasma. 

Expression ([6]) implies that if the monopoles affect the 
gluonic condensate then the monopoles should also make 
influence on the trace anomaly and, as a consequence, on 
the equation of state of the gluon plasma. Translated to 
the lattice language p4)) . this statement means that the 
thermal monopoles contribute to equation of state if they 
affect the expectation value of the plaquette action Sp. 

At zero temperature the correlation of the monopoles 
with the action density was indeed established in 
Refs. [l^ll^llO]- Although this result was obtained in a 
cold vacuum, it provides us with a hint to suspect that 
the monopoles may contribute to the equation of state at 
finite temperature as well. 

In the continuum limit, the averaged value of the action 
density at the position of monopoles can be calculated 
with the help of the normalized correlator between the 
monopole trajectory j^{x) and the non-Abelian action 
density, 



(1^) • {{Sp)t (5p)o) . (14) (5)„on(x, x') 



2g2 ^"^^ G f_ii/)n\on{x, X ) 



(16) 
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The points x and x' are very close to each other. As we 
will see below the lattice provides a natural regularization 
of the correlation function in Eq. (fT6|) . 

The contribution of the monopoles (fT6|) to the gluon 
energy density (or, equivalently, to the gluonic conden- 
sate) can be subdivided into the electric and magnetic 
parts, respectively: 



{S)mon — {Sm)i 



mon ■ 



(17) 



The magnetic part of the action density can naturally be 
expressed as follows [28| : 



{SM)monix, x') = 



1 



{Jlix)) \2g2 



/ 7r^Trb-^(a;)G^.(x')]'Vl8) 



where 



G^v{x) 



TipGapix) . 



(19) 



If the monopole is static, ~ (5^4 then only spacelike 
(magnetic) components of the field strength tensor Gy 
with i,j = 1, 2,3 contribute to the correlator (fT5|) . 

The excess of the chromoelectric action around the 
monopoles with respect to the chromoelectric vacuum 
expectation value in the bulk can be evaluated with the 
help of the following normalized correlator [2^ : 



{Si 



i{x,x') 



1 



{jI{x)) 



iTr[j^(a;)G^.(x')]' 



.(20) 



If the monopole is static, ~ (5^4 then only timelike 
(electric) components of the field strength tensor G4i 
with i ~ 1,2,3 contribute to the correlator (|20l) . 



IV. LATTICE RESULTS 
A. Lattice observables 

The continuum definitions and (PO)) can be trans- 
lated to the lattice as follows. First, we locate the 
wrapped monopole trajectories using a standard pre- 
scription described in Appendix |^ The monopole cur- 
rent is defined on the links {a;,/i}* of the dual lattice 
and takes integer values, jfi{x) = 0,±1, . . .. Usually the 
multiple-charged currents with |j;^(a;)| ^ 2 are very rare, 
so that one can set |j^(a;)| — jf^{x) with a high accuracy. 
Each link {x, fi}* of the dual lattice corresponds to the 
3-cube Cp(x) of the original lattice. Then Eq. (fT8)l in 
continuum has the following lattice analogue [28[: 



(5 



M 



(21) 



X fl—l 



y X /J, — 1 



Pedc^ix) 



where the sum over plaquettes is taken over all six faces 
P of the 3-cube Cf_i{x). 

The lattice definition (|2ip corresponds to a minimal 
splitting between the points x and x' in the continuum 
expression p8|) . The center of the monopole coincides 
with the center of the 3-cube Cp,{x) which surrounds the 
monopole. On the other hand the monopole contribu- 
tion (j2ip to the magnetic part gluon action density is 
defined by the average of the plaquettes at the faces of 
the cube. Thus, there is a minimal distance between the 
monopole position and the point where the gluonic action 
density is measured. Thus, the lattice formulation (PT|) 
of the magnetic correlator (jlSp provides us with a natu- 
ral definition of the splitting between the points x and x' . 
The minimal distance equals to the half lattice spacing a: 



(22) 



The lattice analogue of the chromoelectric correla- 
tor pH) is j2a]: 



{SE)r 



EE^"(-) 

y X /J,— 1 



EE^-^(-) 
- E 



(23) 



24 ^ 1^-2^^^- 

Pedv[c^(x)] 



1. 



where one of the sums is taken over 24 plaquettes P which 
form a manifold V[Cf_i{x)] corresponding to the 3-cube 
Cf^{x). This manifold is represented by a set of all pla- 
quettes P which satisfy the following two conditions: 

(i) each plaquette have one, and only one, common 
link with the cube C^{x); 

(ii) each plaquette is lying in the planes, defined by the 
vectors /i (the direction of the monopole current) 
and p (the direction of the link defined in the pre- 
vious condition). 

One can estimate a minimal distance between the glu- 
onic operator Sp and the position of the monopole in the 
lattice formulation (|23|) of the electric correlator ((20|) : 



V2 



(24) 



In order to make this estimation we assumed that the 
monopole trajectory is static. This assumption is jus- 
tified because as the temperature increases the zeroth 
Matsubara component for all bosonic fields - including 
the monopole trajectories - becomes dominant. 

The excess of the magnetic and electric parts of the 
gluonic action due to the presence of the monopoles is 
given by the following formulae, respectively. 



SSa = {Sa), 



{S, 



A = M, E . 



(25) 
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where (5'A)vac niagnetic (for A = M) and electric 

(for A = E) action densities in the vacuum: 




Here the sums go over spatial (Pspat) and temporal 
(-Ptemp) plaquettes, respectively. 



A ^ 









/ 


\^ 

.01 




1 

0. 






f r 




02 






1 



FIG. 1: The excess of the magnetic component &Sm of the 
gluonic action (|25p as the function of the distance r, Eq. (|22p 
from the monopole and the temperature. All quantities are 
given in units of the critical temperature Tc. 



B. Observation of electromagnetic duality 

In Fig. [T] we plot the excess of the magnetic compo- 
nent 8Sm of the gluonic action l[25|) over the correspond- 
ing vacuum expectation value: results are reported in 
dimensionless units and as a function of the distance r 
from the monopole center and of the temperature. We 
include into this figure all available data sets which are 
described in Appendix A in more details. Notice that 
available data refer to a single value of r for each fixed 
value of the UV cutoff, so that results reported in Fig. [T] 
have been obtained by collecting together data from dif- 
ferent simulations performed at different values of the lat- 
tice spacing. Therefore the dependence of 8Sm on r, as 
reported in Fig.[l] cannot be disentangled from a possible 
dependence of 8Sm on the UV cutoff: more simulations 
in which the non-abelian action density is measured at 
more than one value of r for each lattice spacing would be 
necessary to completely address this issue. For that rea- 
son, our present discussion about the dependence of (5S'j\/ 
on r should be considered as preliminary and qualitative. 

The excess of the magnetic action near thermal 
monopoles turns out to be a positive quantity at all 
distances r from the monopole and at all tempera- 
tures, moreover it decreases as a function of r. This is an 
unexpected result because at zero temperature the con- 
densed monopoles have an excess density of magnetic ac- 
tion which, even if still positive, increases as a function of 



r, i.e. it is minimum at the monopole center [30|. As the 
monopoles from the condensate evaporate into the ther- 
mal monopoles, one could expect that in the deconfine- 
ment phase the excess of the action near monopoles has 
the same qualitative features as below T^- This expecta- 
tion turns out to be not the case, however. The depen- 
dence on T is much milder than that on r: the excess of 
the magnetic action is almost temperature- independent. 
The same qualitative features are also seen for the electric 
part of the action, which we show in Fig. [21 




FIG. 2: The same as in Fig. [T] but for the electric part 5Se- 
The distance is defined by Eq. ([24]). 

It is interesting to note that the electric and magnetic 
contributions practically coincide with each other, 

(£;2) « (s2^ ^ (28) 

where 

E' = \gUx) and B' = )gI{x), (29) 

are the electric and action densities, respectively [in our 
units Eq. ([^5)1 means bSu ~ &Se\- On the contrary, for 
a purely "chromoelectric" object, like quark, one expects 

(£;')qua,k^O, and (i?')q,ark'-0, (30) 

while for a purely "chromomagnetic" object, like a 
monopole, one should get an opposite relation, 

(i?'),no„~0, and (i?'),„o„y6 0. (31) 

In Fig. [3] we graphically demonstrate the validity of 
Eq. ([25|) by plotting the overlap of the Figs. [U and [1 
For convenience, we show the both contributions to the 
energy density multiplied by the distance r. Basically, 
Fig. [3] represents a scaled projection of Figs. [1] and [2| 
onto the distance-action plane. This figure illustrates 
the amazing fact that the electric and magnetic actions 
are the same around the "monopoles". Indeed, the fact, 
that the electric and magnetic contributions are almost 
the same is unusual because from the first principles we 
would expect that the magnetic component of the glu- 
onic action near the monopole is much stronger than the 
electric part of the gluonic action. 

Even if the equality of electric and magnetic contribu- 
tion is observed only on average quantities, it is tempting 
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to make the hypothesis that thermal abehan monopoles 
are associated to non-abehan field configurations possess- 
ing EM duality, i.e. selfdual objects. Numerically, the 
deviations from the EM duality do not exceed 3.5%. 

The objects which carry equal amounts of electric and 
magnetic components are non-Abelian dyons. The (clas- 
sical) fields of dyons satisfy the selfduality relation 



G tiu 



1 



(32) 



so that electric and magnetic action densities ((29|) of the 
dyons are the same. 



E. 



clyoii 



dyon • 



(33) 



This equation is in line with our numerical observa- 
tion (pg]) . 

The simplest selfdual configuration is a (static) 
Bogomolny-Prasad-Sommerfield (BPS) solution to the 
classical equations of motion of SU{2) Yang-Mills the- 
ory [IH: 



g r 

1 

At = -Hr)-, 
g r 



fir) 



1 



1 



sinh ) 



, (34) 



h{r) = -(rcothr - l) .(35) 



Here the distance coordinates r = |a;| and x are scaled by 
an arbitrary factor, r —^ ■ r, to make a dimensionless 
quantity. 
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FIG. 3: Illustration of the dyonic nature of the objects: The 
magnetic (the filled circles) and electric (the filled squares) 
contributions to (the square root of) the action density mul- 
tiplied by the distance r as the function of r. The statistical 
error bars are much smaller than the symbols used. 

In order to characterize both the singularity at the 
origin of the thermal monopole (dyon) and the observed 
EM duality we fit the magnetic and electric excesses of 
the action density by the function 



r 



(36) 



We performed three different fits. We fitted the elec- 
tric action density, £ = E and, separately, the magnetic 



action density, i = M, and finally, we made the fit of 
the both simultaneously, i = EM. In all these fits we 
have averaged our data points corresponding to different 
temperatures T and the same distances r (this is a legiti- 
mate operation since the data are almost T-independent 
as one can see from Fig. [3]). The best fit parameters are, 
respectively: 



Ce = 0.243(2), 
Cm = 0.247(2), 
Gem = 0.246(2). 



(37) 



These values are very close to each other, so that the du- 
ality works with a high precision. In Fig. H] we plot 
the joint Rt £ — EM in order to illustrate the quality of 
the selfduality relation Thus, in this paper we refer 
to these objects as both monopoles and selfdual dyons 
simultaneously, because they are identified as monopoles 
by their magnetic charge, while their non-Abelian con- 
tent shows signatures of the dyonic selfduality. 
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FIG. 4: The fit of the excesses in the magnetic (the blue 
circles) and electric (the red squares) action densities near the 
monopole as the function of the distance r from the monopole 
(shown both in units of T^^ and in Fermi). The error bars 
and the magnitude of the temperature dependence is much 
smaller than the size of the symbols. The data are fitted by 
the function ((36]). The fit is shown by the dashed line. 



Note that (smooth) dyonic configurations were found 
in Ref. [l^l using a method of cooling of the gluon con- 
figurations and/or utilizing the low- lying fermion modes 
to filter the positions of the dyons. In this paper, assum- 
ing that thermal monopoles are indeed associated with 
dyons, we report the direct observation of dyons in dy- 
namical (uncooled) configurations. In fact, the field of 
the classical BPS dyon (|35p is a smooth function at 
the origin, r = 0. On the contrary, our data suggest that 
the action of the thermal dyons behave as singular as 
l/r'*. One would conclude that quantum thermal dyons 
are not resembling the smooth classical BPS solutions at 
all. On the other hand, as already stressed above, our 
measurements have been taken at the scale of the lat- 
tice spacing, so that the dependence on r cannot really 
be disentangled from that on the UV cutoff and our dis- 
cussion about the 1/r^ behavior should be regarded as 
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only qualitative: more refined measurements should be 
performed in order to clarify the issue. 

C. Contribution of thermal monopoles to pressure 

Finally, let us discuss the contribution of thermal 
monopoles (or, dyons) to the thermodynamics of the 
plasma. The positive value of the excess of both electric 
and magnetic action densities near the monopoles means 
that their contributions to the trace anomaly 0, Eqs. ([6]), 
(fH)) and, consequently, to the pressure p and to the en- 
ergy density e, Eq. ([8]), are negative quantities [s^ : 

5e^.on{T) < 0, 

^Pmon(T)<0, (38) 
&mo„(T)<0, 

because the /3-function /3(g) in Eq. ([S]) is a negative- 
valued function. We cannot calculate numerical values 
of the monopole contribution to these thermodynamic 
quantities because we should also take into account the 
entropy of the monopoles. This task cannot be completed 
using our method. We can only state the qualitative fea- 
tures of the effect of the monopoles 

V. CONCLUSION 

In this paper we have studied the correlation of ther- 
mal abelian monopoles detected after MAG projection 
with the gauge invariant nonabelian action. The study 
has been performed for QCD with two colors and in the 
quenched limit. We have found an equal excess of electric 
and magnetic action around thermal monopoles, suggest- 
ing that the underlying non-abelian gauge configurations 
- associated to thermal monopoles - may correspond to 
selfdual dyons. The electric-magnetic equality holds with 
a good numerical precision (less than 4%). 

We notice that the presence of the selfdual dyons in 
the thermal plasma is a nontrivial fact. Indeed, one can 
argue [l^ that the thermal medium is populated with 
gluonic monopoles, which become dyons via the Witten 
effect dH as these monopoles pass through self-dual re- 
gions (say, instantons) in space-time. In this case the 
monopole would get a space-dependent electric charge 
which, however, will be vanishing far from the instan- 
tons. In this scenario the electric action density should 
be small compared to the rival magnetic component. On 
the contrary, we have observed exact selfduality indicat- 
ing that - at least at short distances - the electric com- 
ponent is as strong as the magnetic one. Note that the 
electric-magnetic duality may be a property peculiar to 
thermal monopoles observed in the deconfined phase of 
Yang-Mills theories. 

We have found that thermal monopoles (dyons) make 
a negative contribution to the pressure and to the energy 



density of the gluon plasma. 

Most importantly, certain signatures of the dyons may 
be observed experimentally in the heavy-ion collisions via 
the so-called chiral magnetic effect. Indeed, since the 
dyons are characterized by nonzero topological charge 
density, they are playing a role of "seeds" of topologi- 
cal charge fluctuations in the plasma. In a strong mag- 
netic field the topological structures in the gluonic vac- 
uum leave their footprints via certain CP-odd correla- 
tions like, for example, the generation of the electric cur- 
rent due to the chiral magnetic effect [2l[. Preliminary 
signatures of such a current were already observed exper- 
imentally by the STAR collaboration at RHIC ^2^] (this 
effect was also observed in the lattice gauge theory in 
Ref. [13). 

In order to clarify further the nature of these objects 
we are planning to calculate correlations of the thermal 
monopoles with the local topological charge density. We 
are also going to measure cross-correlations of fluctua- 
tions of electric and magnetic action densities around the 
thermal monopoles, and extend our study of correlation 
functions to larger distances from the monopoles. 



Acknowledgments 

This work has been started during the workshop "Non- 
Perturbative Methods in Strongly Coupled Gauge Theo- 
ries" at the Galileo Galilei Institute (GGI) in Florence 
and was supported by the STINT Institutional grant 
IG2004-2 025, by the RFBR 08-02-00661-a, DFG-RFBR 
436 RUS, by a grant for scientific schools NSh-679.2008.2, 
by the Federal Program of the Russian Ministry of Indus- 
try, Science and Technology No. 40.052.1.1.1112 and by 
the Russian Federal Agency for Nuclear Power. Numeri- 
cal simulations have been performed on a computer farm 
in Genova provided by INFN. 



APPENDIX A: NUMERICAL SIMULATIONS 

Our numerical simulations were performed in SU (2) 
lattice gauge theory with the standard plaquette action. 
The basic setup is the same as the one adopted in Ref. fJ] . 
We used various lattice geometries Li^ x Lt with varying 
temporal extension of the lattice Li = 4, . . . 13, and fixed 
spatial extension Ls = 64. We performed simulations 
at different values of the lattice gauge coupling (3 — 
2.7, 2.86, 2.93, 3.00, 3.05. 

The physical value of the lattice spacing was deter- 
mined from the relation 

a(/3)Ai = R{P)m , (Al) 

where R is the two-loop perturbative /3-function, while 
A is a non-perturbative correction factor computed and 
reported in Ref. [l^l- The mass scale A^ is related to 
the critical temperature [111, Tc = 21.45(14)Ai,, and 
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to the zero-temperature string tension a via Tc/ ^/a = 
0.69(2) The set of our values of (5 corresponds to five 
lattice spacings in the interval a = (0.015 ... 0.045) fm. 
The set of temporal extensions Lt and the lattice spacings 
a determines a grid of temperature values T — l/{Lta) 
which is rather wide, T/T^ = 1.095, • • • , 5.637. 

The currents of the Abclian monopole were identified 
by the standard De Grand- Toussaint construction [3^ 
after performing an Abelian projection in the Maximal 
Abelian gauge. This gauge is defined by the condition of 
the maximization of the following functional with respect 
to gauge transformations: 

Y^Retr[U^{x)a^Ul{x)az] (A2) 

tJi,X 



The maximization is achieved by an iterative combina- 
tion of local maximization and overrelaxation |37i | . 
In the Maximal Abelian gauge the monopole currents 
are defined as location of sources of the magnetic fields 
in the diagonal components of the gauge fields [1^: 

ZTT 

Here O^i, is a compactified plaquette angle, — tt < 9^,, ^ 
TT, defined as follows: 

dt„y = dt,u + 27r7i^^ , n^^ e N . (A4) 
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